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Abstract 

We show how studying leading singularities of Feynman diagrams, when all momenta are com- 
plex, gives a simple way of writing multi-loop and multi-particle scattering amplitudes in M = 4 
super Yang-Mills. The simplicity of the method is equivalent to that of the quadruple cut technique 



introduced in hep-th/0412103 at one-loop. The new technique only involves the computation of 
residues and the solution of linear equations. In our technique both parity even and parity odd 
pieces of a coefficient are computed simultaneously and it is only at the end that a separation can 
be made if desired. We explain the procedure via examples. The main example, which we compute 



in detail, is the five-particle two- loop amplitude first given in hep-th/0604074. Another feature 
of our method is that the helicity structure of the amplitude only enters in the inhomogeneous 
part of the linear equations. In other words, the homogeneous part is universal. We illustrate this 
feature by presenting the linear equations which determine a large class of terms for MHV and 
next-to-MHV six-particle two-loop amplitudes. 
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I. INTRODUCTION 



The idea of determining the S'-matrix of a theory from the structure of its singularities 
was intensively studied in the 60's The main idea was to consider the 

S'-matrix as an analytic function of the kinematical invariants with singularities determined 
by physical input. Most efforts were directed towards the study of hadronic interactions. 
Nowadays we know that understanding analytically the strong interactions is a very hard 
problem. 

A more modern version, similar in spirit, was developed in the 90's, where extensive use 
of branch cut singularities was shown to be a powerful tool in the computation of scattering 
amplitudes in Yang-Mills theories. These techniques came to be known as the unitarity 



based method 
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More recently [14J, [15|, H6J , the highest codimension singularities, called leading singu- 
larities already in the 60's, were shown to be a powerful tool in computing amplitudes. 
At one-loop, the problem of computing amplitudes is reduced to that of computing tree 
amplitudes for theories where no triangles or bubbles appear. One such theory is M = 4 
super Yang -Mills (SYM) and it has been hypothesized that so too is Af = 8 supergravity 
2( 
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In this paper we show that at hig 
has only been partially unleashed 15 



rer loop orders the power of the leading singularity 



161 ]. By carefully studying all leading singularities we 



propose a method which reduces the computation of multi-loop and multi-particle ampli- 
tudes in Af = 4 SYM to the computation of residues (which end up being related to tree 
amplitudes) and to the solution of linear equations. 

The discontinuity across the leading singularity at one loop is computed by collecting all 
Feynman diagrams that share four given propagators and then cutting the propagators to 
turn the sum into simply the product of four tree-level amplitudes. Cutting means removing 
the principal value part of a Feynman propagator, i.e. l/(i 2 + ie) — > 5^(£ 2 ). Thus the 
integral over the four-vector I is completely localized by the four delta functions. The value 
of the integral for each solution is given by the jacobian of the change of variables from I 
to the argument of the delta functions evaluated at £*. The final answer is then the sum over 
of all contributions. In massless theories there are two solutions I* and i\ . In general, 
the product of the four tree-level amplitudes gives different answers at the different points 
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{£*,£*}■ An important exception is when the number of particles is four. In this case, the 
product of amplitudes is equal to A tTee st for both values of £. 

In general, the support of the delta functions is outside the region where I is a real vector. 
This means that the computation of the leading singularity 1 is more naturally interpreted 
in terms of a contour integral in C 4 where i is now a complex vector. There are two distinct 

n 

leading singularities and the prescription of [14j is equivalent to choosing a contour that picks 
up both residues. A natural question is what the role of each isolated leading singularity is. 
Naively, one could try to expand the one-loop amplitude in scalar boxes and then compute 
their coefficients by matching the residues at a given leading singularity. If one did that one 
would find different answers for the same coefficient which would be a contradiction. 

In this paper we show that the resolution to this puzzle is that at one-loop, reduction 
formulas that express e.g. pentagons in terms of boxes, derived in 
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231 ] . are not valid 



for generic complex integration contours. This means that in order to write an expression 
in terms of scalar integrals which reproduces all singularities of Feynman diagrams correctly 
one has to allow for higher point integrals. This explains why four particles is special and 
why the first non-trivial case is that of five particles, where one should allow for a pentagon. 



At higher loops, the original approach of 15] is to sum over all solutions. In [16], it was 



stated that one should be able to work with individual singularities, but this was in the 
context of four-particles where there is no distinction. 

In this paper we show that requiring agreement between Feynman diagrams and scalar 
(or generalized scalar) integrals on each individual leading singularity provides enough linear 
equations to determine higher-loop and multi-particle amplitudes. 

In an nutshell, one starts collecting all Feynman diagrams with chosen propagators that 
give rise to topologies with only boxes. By writing an ansatz of scalar integrals (i.e., with 
numerator one) one requires agreement at all leading singularities independently. This gives 
a set of linear equations. Sometimes the system of equations does not have solutions which 
means that generalized scalar integrals (i.e., with numerators that are propagators with 
negative power) must be added. The numerators serve as zeroes to cancel poles so that the 
new integrals have zero residue on leading singularities where they are not needed. 



1 Here and throughout the paper we will abuse terminology and refer to the discontinuity across a given 
leading singularity as the leading singularity itself. 
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We explain the method via examples. Our main example is the five-particle two-loop 



amplitude. This amplitude was first computed in 
requires solving linear equations in two variables! 



24j . In this case, our computation only 



An interesting new feature already found in [25j, as compared to the four-particle am- 
plitude, is that after normalizing by the corresponding tree-level amplitude one finds terms 
that are parity even and some that are parity odd. Using our technique, it becomes very 
transparent why the answer does not have definite parity. The reason is that leading sin- 
gularities come in pairs. For real external momenta, one of them would be located at the 
complex conjugate value of the other. For four particles, the value of Feynman diagrams on 
each leading singularity is the same and hence one gets a parity even answer. For five or 
more particles, the values differ and one gets complex solutions. 

It is important to mention that the actual helicity structure of the amplitude only enters 
as the inhomogeneous part of the linear equations that determine the coefficients of the 
integrals. In this sense, the homogeneous part of the linear system of equations is universal. 



Since for five particles, all non-trivial helicity configurations are MHV (or MHV) the feature 
just described is not very surprising. This is why we present the linear equations which 
determine a large class of terms contributing to six-particle MHV and next-to-MHV two- 
loop amplitudes. Solving the equations could allow a comparison to the recent computation 



27|]. 



of the parity even part of the six-particle MHV amplitude in 

This paper is organized as follows. In section II we establish our conventions. In section 
III we give a detailed explanation of the leading singularity technique at one-loop. The 
discussion of how the scalar basis must be extended is given for five particles since the results 
there can be directly applied to the two-loop case. In section IV we introduce the leading 
singularity technique at two loops in the simplest case, i.e., the four-particle amplitude. In 
section V we present the main example of the paper, the five-particle two-loop amplitude, 
in complete detail. In section VI we present the linear equations for coefficients in MHV 
and next-to-MHV six-particle two-loop amplitudes. In section VII we give conclusions and 
future directions. 
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II. PRELIMINARIES 



Scattering amplitudes of on-shell particles in Af = 4 SYM with gauge group U(N) can be 



written as a sum over color-stripped partial amplitudes using the color decomposition 

y. 



id, 



, 3(]] . Each partial amplitude admits a large N expansion. More explicitly, 

A(l,2,...,n) = 5^( Pl +p 2 + ...+ Pn ) Tr (T ai T a2 ...T a ") A n (l, 2, . . . , n) 

+ permutations + . . . (1) 

where the sum is over non-cyclic permutations of the external states (cyclic ones being 
symmetries of the trace) and the ellipsis represents terms with double and higher numbers 
of traces. A n may be expanded in perturbation theory and we denote the L-loop planar 
partial amplitude by An . We also use A^ 1 = A^ ee . 
Our conventions are: 

• A tree-level MHV amplitude has two particles of negative helicity and the rest of 
positive helicity. 

• A null vector p^ is written as a bispinor as p ail = X a Xa- 

• The Lorentz invariant inner product of two null vectors p^ and is given by 2p ■ q = 
(A p , Xg)[X p , Xg). We also use (X p , X q ) = (p, q) and [X p , X q ) = [p, q]. 

• All external momenta in an amplitude are outgoing and usually denoted by hi. 

• Some useful Lorentz invariant combinations are Sij — (ki -\- kj^ } tiji — (fcj -|- kj -\- fc/)^, 
and [a\b + c\d) = [a, b](b, d) + [a, c](c, d). 

III. LEADING SINGULARITY AT ONE-LOOP 

Any one-loop amplitude in Af = 4 SYM may be expressed both as a sum over Feynman 
diagrams, and also in terms of scalar box integrals Q]: 

A^ = ^ {1-loop Feynman diagrams} = ^ B i x A K l , K 2, K s, K l) ( 2 ) 

x 

where the second sum is over all partitions X of {1,2, ... ,n} into four non-empty sets, Kf 
equals the sum of the momenta in the i th subset of partition I and the Bj are coefficients to 




m-1 



= B x 



K 



FIG. 1: Computation of a coefficient using the leading singularity of a box. The lines circling the 
propagators represent the T 4 contour of integration. The left hand side of the figure represents the 
sum of all 1-loop Feynman diagrams - note that only those Feynman diagrams that contain the 
displayed propagators actually contribute to this particular contour integral. 

be determined. Since we are working with color ordered Feynman diagrams, one considers 
only partitions that respect the color ordering. 
Scalar box integrals are of the form 

dH 1 



I{K X ,K 2 ,K Z ,K A ) :-- 



(2tt) 4 i\i - K x )\t -K x - K 2 ) 2 (£ + K 4 ) 2 



(3) 



Since scalar integrals are known explicitly (see e.g. 28(] ) , the problem of computing any 
one-loop amplitude is reduced to that of determining the coefficients B%. 

The amplitude, denned in terms of Feynman diagrams, possesses many singularities. 
Some of them are branch cut singularities and comparing them on both sides is a way of 
obtaining information about the coefficients and in many cases it allows their determination 

Q 

like for MHV amplitudes [7|j. The main difficulty is that a given branch cut is generically 
shared by several scalar integrals and therefore several coefficients appear at the same time 
and the information must be disentangled. 

It turns out that the highest codimension singularities of Feynman diagrams, called lead- 
ing singularities, receive contributions from a single scalar box integral. Thus, using leading 
singularities of Feynman diagrams one can determine all coefficients, Bj, one by one fl^ |. 

If we let fi(£) with i — 1, ... ,4 correspond to the four factors in the denominator of (J3J), 
then the leading singularity is computed by replacing 1/ fi(£) by 6(fi(£)). Applying this to 
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both sides of (EI) one finds that 



E det 



dfi 



-1 



Multiplet i=l 



5^ x det 



dfi 

din 



(4) 



where the sum over means a sum over all solutions to the equations fi(£) = for i = 
1, . . . , 4. In general there are two solutions. The second sum on the l.h.s. is over all members 
of the M = 4 supermultiplet as choices of internal particles. 

When the number of external particles is four, i.e. n = 4, both solutions give a jacobian 
factor equal to l/{st) and the sum over the product of tree-level amplitudes equals A tTee st. 
This implies that B = A tree st. For five or more particles the sum over tree amplitudes is not 
the same for both solutions. In fact, for five particles the product of tree amplitudes vanishes 
in one solution and it gives A tree st in the other. Here s = (Ki + K 2 ) 2 and t = (K 2 + K 3 ) 2 
with KiS as in ([3]). Using this in (j4j) one finds that B = A tvee st/2 which is the correct 
answer 30j. Note the factor of half coming from the fact that on the r.h.s. of (jlj) one gets 
a factor of two. 



A. Sharpening The Leading Singularity 

The solutions, £*, to the equations fi(£) = are complex in general. This makes the 
mathematical interpretation of the leading singularity more transparent if defined as a con- 
tour integral in C 4 . The contour integral is obtained by simply taking ([3]) and allowing £ to 
be a complex vector. The contour of integration has the topology of a four-torus, T 4 = (5' 1 ) 4 , 
around each isolated singularity, and it is given by T = {£ : \fi(£)\ = e{\ with q some small 
positive number. The integral is computed by a generalization of Cauchy's theorem to higher 



dimensions (see for example section 5.1 of 3J|). This can be seen by a change of variables 



to local coordinates Zj's in terms of which the singular point is located at the origin, i.e. 
z % — fi{£) — 0. Now it is clear that the contour of integration, T = {£ : \zi\ = e^}, when 
projected on each ,2 r plane is just a circle of radius e^. The integral becomes 2 



Here and throughout the paper we use the convention that the symbol § contains a factor of l/(27ri). 
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A natural questions that comes from this formalism is how to compare the sum over 
Feynman diagrams and the sum over scalar integrals when evaluated on each leading sin- 
gularity independently. In other words, we would like to make sense of (jlj) term by term in 
the sum over solutions. The reader might anticipate a problem from the fact that the sum 
over tree amplitudes gives different answers for n > 4. Before giving the resolution to the 
puzzle explicitly for n = 5 let us discuss in detail the four-particle case. 

1. Four Particles 

For four particles there is a single planar configuration with no triangles or bubbles. This 
is depicted in figure [2]A The location of the two leading singularities is: i a a = aAo and 
4d = aAi 4) Ai 1} with a = [l,2]/[4,2] and a = (l,2)/(4,2). 

It is straightforward to compute the product of the four three-particle tree-level ampli- 
tudes on each of the solutions. In each case we find the same answer A tree st. In order 
to reproduce this with scalar integrals we start with a single box and a coefficient to be 
determined. 

Comparing both sides evaluated on each of the two leading singularities we find the same 
condition: 

B = A tree st. (6) 

This indeed gives the correct answer. As we will see, the fact that we got the same condition 
by requiring the two sides to match on each leading singularity independently is only a 
peculiarity of n = 4 and it will not be true for any n > 4. 

2. Five Particles 

Five particles is the first non-trivial example where the conditions coming from the two 
leading singularities are different. At this point there seem to be a puzzle since requiring 
agreement on each singularity independently leads to distinct coefficients for the same scalar 
box integral which is indeed a contradiction. 

The resolution to this puzzle is simple but subtle. It has to do with the fact that when 
momenta are complex, the integrand of loop amplitudes is invariant under the complexi- 
fication of the (double cover) of the Lorentz group, i.e. SL(2,C) x SL(2, C). However, 
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(A) (B) 

FIG. 2: Sum over Feynman diagrams with no triangles or bubbles. (A) Unique four-particle 
configuration. (B) Unique five-particle topology with some choice of external labels. For this 
particular choice of external helicities there is a single configuration of internal helicities. 

the integral might be invariant under different subgroups of SX(2, C) x SX(2, C) depending 
on the choice of contour 3 . The physical choice of contour requires the loop momentum to 
be a real vector and hence it breaks the group SL(2, C) x SL(2, C) down to the diagonal 
SL(2, C)d which is a double cover of the physical Lorentz group SO(3, 1). In particular, on 
the physical contour, scalar integrals are parity invariant. This is what fails on the individual 
T 4 contours. 

Let us illustrate this by computing a five-particle MHV amplitude at one-loop. 

Once again there is a single topology of a planar configuration of Feynman diagrams with 
only boxes. This is shown in figure 03 with a particular choice of labels. In this case, the 
location of the two leading singularities is 

<£> = A.« (a! 1 ' + ||a«) , ,2 = (a« + Ma?) X« (7) 

Evaluating the product of the four tree-level amplitudes is again straightforward and 
gives 

4 4 

_ „ ' .-, _ _. z n °- ^ 

, n Multinlet i = 1 ^oi 

q=q\ 



Multiplet i=l 



q=q 



(l) Multiplet i=l 

We want to reproduce this behavior using scalar integrals. The natural candidate is the 
scalar box integral in figure [3]A.. It is easy to see that on each of the leading singularities 
the scalar box integral gives the same residue, i.e., 1/(512523). If this box were the only 
contribution on the scalar integral side we would find that by comparing to the first equation 



A similar situation was encountered in the derivation of MHV diagrams from twistor string theory in 32 1. 
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(A) 




FIG. 3: (A) Scalar box integral for five particles. The particular choice of labels corresponds to 
the integral /( a ) :1 in the text. (B) Scalar pentagon integral. This integral is denoted by . As in 
the rest of the paper, external momenta are taken to be outgoing. 

in (jHJ) the coefficient would be B = v45 ree si2S23, which is twice the known answer. If, instead, 
we use the second equation in (jHJ) we would find B = which is also wrong. 

This is the puzzle mentioned earlier. In order to discover the resolution, let us assume we 
did not know that M = 4 SYM amplitudes can be written purely in terms of scalar boxes 4 . 
The natural starting point to reproduce the behavior of the collection of the Feynman 
diagrams is the scalar box integral in figure [31 However, as we have seen, this is not enough. 
Therefore we need to expand the basis. The only other scalar integral that can contribute 
to the same leading singularities, but with different residues, is a pentagon (see figure [3]B)- 

Denoting the coefficient of the pentagon by C we find that reproducing (jHJ) means 

B + W ^ ? =A^ Sl2 s 2s , B + w ^ = . (9) 

where (q + k 5 ) 2 is the only propagator of the pentagon not shared by the box. 

These equations can easily be solved. A convenient way to express the solution is obtained 
by making the following definitions 

,_ f, , {i + 2,i + 3)[i + 2,i] y 1 ~ _ ( (z + 2,z)[i + 2,z + 3] y 1 
Pi '~ V (i + l,i + 3)[i + l,i\J ' Pi '~ V + (i + l,i)[i + l,i + 3]J ' 1 ' 

Note that for real momenta, is the complex conjugate of This will be important when 
identifying odd and even pieces under parity. 



In dimensional rcgularization this is true only to order e° except for n = 4 when is true to all orders in e. 
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The solution to the equations in flU]) is 



,i ree Sl2S23/35 n /1 tree S 5lSl2S23 



g _ _^tree^J^_££^o q _ ^tree "Qi"i^"^ 

The attentive reader might anticipate another possible contradiction: If we have fixed 
the coefficient of the pentagon by using the leading singularity in figure [2B, it is hard to 
believe that it will simultaneously solve the equations coming from figure [2]B after a cyclic 
permutation of the labels 5 . For this to be true the same coefficient of the pentagon should 
work in all e., C/A tTee in (JTTj) must be invariant under cyclic permutations of the 

labels of the external particles! 

Indeed, an explicit computation reveals that 

q C _ s i,i+l s i+l,i+2Sj+2,i+3 

-A*"~ Pi- Pi 

is the same quantity for all i £ {1, 2, 3, 4, 5}. How to write C in a manifestly invariant form 
will become clear in the next section. 

We are ready to write down the final answer for the amplitude. In order to compare with 
the known result it is convenient to separate contributions into parity even and parity odd 
pieces. Note that the coefficient of the pentagon has definite parity; it is parity odd. The 
coefficient of the box does not have definite parity. In order to decompose it we write in the 
numerator of B, /3 S = (/? 5 — /?s)/2 + (/3 5 + /3 5 ) /2. The amplitude is then given as 

A (1)MHV _^ i{ajA J / ^ + ^ Sl2S23 /H:l ] + C /«. (13) 



where is a scalar box integral with three massless legs given by the momenta of the i th , 
(i + l) th and (i + 2) th particles while 1^ is a pentagon with all massless legs as shown in 
figure ©B. 

We claim that this is the correct answer. In order to prove it we have to go back 
to the usual contour of integration where I is a real vector and where integrals must be 
dimensionally regulated. In that case we have to show that (JTBl reduces to 

,(1)MHV 



^treeMHV 



Dim. Reg. cyclic 



E (5* • (14) 



Note that the pentagon is invariant under a cyclic permutation of the labels. 
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Comparing our answer (|T3|) with the answer for real I and in dimensional regularization 
we conclude that the scalar box and pentagon integrals when dimensionally regulated must 
satisfy the following identity 

j(6) = -^2 /(«)*. (15) 



In other words, a pentagon can be written as a sum of five boxes with 



particu 



ar coefficients. 



This is an example of what is known as a reduction formula [211. 221 . |23| . Indeed, one 
can check that the coefficients in ( 115i) agree with those obtained in 23[ by using differential 
equations, i.e., 

— + — = {oci-2 - oti-i + oii- a i+ i + ai + 2)aia 2 a3a4a5Si+i,i+2Sj + 2,j+3 (16) 

s i,i+l 

with ai = Sj + i i i + 2Sj + 2,i+3/H and S = \/ — s i2 s 23 s 34 s 45 s 5i- If we had known the existence of 
reductions formulas but not their form then the leading singularities would have given an 
elementary way of finding them 6 . 

Let us conclude this section by explaining why reduction formulas do not hold on the 
T 4 contours used to compute independent leading singularities. First observe that the 
choice of one such contour breaks the Z2 symmetry that exchanges the two factors in 
SL(2, C) x SL(2, C) which is the symmetry of the integrand. This is nothing but a parity 
transformation. The reason this affects the reduction formulas is that their derivation relies 



on the identity 



22 



f P>* R u R p R a 



with T a real contour 7 . 

A simple way to prove this is by noting that 

" = I " w+imiTmrnw = AiK+ MBt + AsR ° (18) 

for some scalar functions Aj. This is obviously true by Lorentz invariance including parity 
invariance. In general any four vector, in particular, 7^ can be expanded in a basis of vectors 
given by i?^, R%, R% and & vpa R\ v R2 P Rj,a- The latter does not contribute in (fl8l) because it 
is not parity invariant. 



In comparing with [23( , we had to shift the index i because what we call I^ 11 is called I^- 1 
7 If any of the i?f vanishes this integral is divergent. The analysis in dimensional regularization was 



performed in 23(. Out T 4 contour renders all integrals finite and this is why we do not need any regulators. 
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Now it is clear that on a given T 4 contour corresponding to a single leading singularity 
the integral is not parity invariant and ffT8l) does not hold due to the presence of the extra 
vector & vpa Ri v R2pR? J(J in the expansion of J M . This also shows why in the original quadruple 



cut technique introduced in 



14j . which is equivalent to summing over both contributions, 



reduction formulas are valid. Summing over both contours preserves parity since a parity 
transformation corresponds to exchanging the two T 4 's. 



3. Higher Point Amplitudes 

Repeating the procedure we applied to the five-particle case to higher point amplitudes 
one should find reduction formulas for higher point scalar integrals since we know that in 
dimensional regularization all n-particle amplitudes are given in terms of boxes. A case of 
special interest in the following is n = 6. We postpone the discussion to section [VI] where 
several coefficient of two-loop six-particles amplitudes are studied. 



IV. LEADING SINGULARITY AT TWO LOOPS 



The physical meaning of leading singularities at higher loops was given in [16]. It was 
shown that Feynman diagrams arrange themselves so that a meaningful T 4i = (S a ) 4L con- 
tour can be defined at L loops. A generalization of what was done at one-loop would require 
interpreting the integral over L loop momenta as a contour integral in C 4i . However, com- 
bining Feynman diagrams that share the topology of a graph with only boxes generically 
gives rise to less than 4L propagators. Therefore defining a T 4L contour is more subtle than 



in the one-loop case. As shown in 



once a partial integration is done on one of the loop 



variables, new propagator-like singularities appear which can be used to define the T con 



tour iteratively. The physical meaning given in [16| is that after the partial integration, one 
produces a set of Feynman diagrams at one loop order less than the original ones and which 
possesses new tree-level factorization channels. Those channels are the hidden singularities 
needed to define the T 4L contour. In this section we restrict our attention to L = 2. We 
start with the n = 4 case in order to clarify the concepts just explained. 



13 



(1) 



(2) 



FIG. 4: Basis of scalar two-loop integrals for four-particle amplitudes in N = 4 SYM. 



A. Four Particles 



We will illustrate the use of the leading singularity technique on a well known case: 
Four-particle two-loop amplitudes. These were first computed in 



25 



331 ] and are given by 



(19) 



where the integrals 1^ and 1^ are shown in figure HI 

We now proceed to reproduce this result using the leading singularities. 
In M = 4 SYM if we sum over Feynman diagrams and consider contour of integrations 
where all legs attached to a one-loop subdiagram are on-shell, such a subdiagram does not 
contain triangles or bubbles. This means that in order to study all leading singularities, we 
only need to consider sums of Feynman diagrams that contain only boxes. In the case at 
hand, these are Feynman diagrams with the topology of a 2-loop ladder. 

A particular choice of external particles is shown in figure Performing the integra- 
tion over the p momentum we find that the product over the four three-particle tree-level 
amplitudes, including the jacobian, gives rise to a four-particle tree-level amplitude |l6 |. 



Therefore, we are left with the sum over one-loop Feynman diagrams shown on the upper 
right of figure [5j Note that we have used that for n = 4 the product of amplitudes gives the 
same answer when evaluated on the two solutions p^ and p^ 2 \ 

The tree-level four particle amplitude has two factorization channels. One of them is in 
the limit when (/C1 + /C2) 2 — > while the other is in the limit when (g — ^i) 2 — > 0. It must now 
be clear how to define the remaining T 4 in order to perform the q integration; one uses the 
original three propagators together with the new propagator l/(g — ki) 2 . On this contour 
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FIG. 5: Computation of the amplitude using the leading singularity of Feynman diagrams. The 
residue of the p integral on any of the two T 4 's represented by the lines circling the propagators is 
a tree-level four particle amplitude. Choosing a T 4 in the q variable which induces a factorization 
in the (q — k\) 2 channel of the tree amplitude one finds the diagram enclosed by dashed lines. The 
residue on the final T 4 is again a four-particle tree-level amplitude. 

the tree amplitude factorizes and gives rise to the diagram in the bottom left of figure 
This is identical to the one-loop case. Once again there are two solutions and q^ 2 \ On 
both solutions the diagrams evaluate to 8 v4 tree . 

On the scalar integral side, we start with an integral of the form shown in figure HI 
After performing the integration over p one finds that the jacobian gives rise to a factor of 
l/(s(q — ki) 2 ). This means that the scalar integral has a non-zero residue on the same T 8 
used for the evaluation of the Feynman diagrams. One finds that the total integration gives 
l/(sH). 

Comparing both sides we conclude that the coefficient of 1^ is s 2 tA i £ ee . 
By using a cyclic permutation of the labels we find that the coefficient of the integral 1^ 
must be st 2 A^ ee . Combining the results we reproduce the known answer ( fl9l) . 

8 The computation also involves the jacobian. In other words, we are computing the full residue. 
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V. TWO-LOOP FIVE-PARTICLE AMPLITUDE 



Five-particle MHV two-loop amplitudes in M = 4 SYM were computed in 24( using the 
unitarity based method. Here we write the known answer and then show how to re-derive 
it using the leading singularity technique. The answer we find comes out in an strikingly 



different form. It is important to make a distinction between MHV and MHV as the two 
amplitudes, when normalized by A\ 
The expression given in [24j is 



l!f ee , are diff eren t. 



A 



(2)MHV 



l^tree,MHV 



+R 



2I^(e)-2s 12 I^(e) + 



S12 
S34S45 



»2A 



i/(a)( 



5+_H 



./(c) ( 



where 



R — ^i234'5l2'523 s 34S45 s 5l/ ^1234) 
$abc — SS1S12 + dSi 2 S 2 3 + &S23S34 — s 45 s 51 + CS34S45, 

ei234 = Aie^pvk^k^kl = tr [^fafafofa] , 



(20) 



(21) 



and 



G 



1234 



det 



\ 



(22) 



/ S12 Sx3 S14 

512 s 2 3 S34 

513 S23 S34 
\ Sl4 s 24 S34 y 

All integrals I^ a \ I {b \ &\ 1^ and 1^ are defined for a particular choice of external 
particles as shown in figure El 



A. Computation Using The Leading Singularity Technique 

We start by following the same steps as in the four-particle case. There are three different 
topologies of diagrams with only boxes. We have to analyze all of them and built an 
expression in terms of scalar integrals which reproduces the behavior of Feynman diagrams 
when evaluated in all possible leading singularities. 

The three topologies are depicted in figure [7] where a particular choice of external labels 
was made. The total set of configurations is then obtained by cyclic permutations of the 
labels. 
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FIG. 6: Basis of two-loop integrals for five-particle amplitudes in N = 4 SYM. Only one choice of 
external labels is shown. The full basis is obtained by considering all cyclic permutation of labels. 

1. First Topology 

Consider the set of Feynman diagrams in figure [7A. This is the first of the three different 
topologies that contain only boxes for five particles. 

Carrying out the integration over a T 4 contour in the p variables we find that for any 
of the two solutions we get a four-particle tree-level amplitude with 1 and 2 as external 
legs (see figured]). In order to continue, we choose a T 4 contour in the q variables that 
induces a factorization of the four-particle amplitude that contains the external particles 1 
and 2. Note that in this case there is a second four-particle amplitude; the one that contains 
external particles 4 and 5. A T 4 which induces a factorization of the second four-particle 
amplitude must also be considered and we postpone its analysis to section IV A 41 Note that 
this second possibility was not present in the four-particle case. 

From the bottom left of figure [S] we find that the problem at hand is identical to the 
one-loop five-particle amplitude discussed in the previous section. The difference between 
the two cases comes in analyzing the scalar integrals. By analogy with the one loop case, we 
start by assuming that the basis contains integrals and 1^ in figure [61 After carrying 
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(B) 



FIG. 7: All inequivalent topologies of sums of Feynman diagrams with only boxes for five particles 
at two loops. A particular choice of external labels is shown. All cyclic permutations of the labels 
must also be considered. 

out the integration over the first T 4 in the scalar integrals, we find that the problem also 
reduces to that of the one-loop case except that the coefficients are multiplied by an extra 
factor of l/si2- In other words, if we denote the coefficients of the two-loop integrals 1^ 
and 1^ by B 2 ~ loop and (7 2 ~ loop respectively, then the equations are 

^2-loop £<2-loop ^2-loop £<2-loop 



S12 si 2 (g (1) + h 



^ tree Sl2S23, 



Sl2 



0. 



(23) 



The solution is obtained by making the substitution B — > B 2 loop /si 2 and C — > C 2 loop /s 

s 12 s 23/35 



12 



in ( TTT1) . i.e, 



^2-loop _ 



tree ' 



^(2-loop 



tree S 5lS i 2 S 2 3 



(24) 



As ~Ps Ps- Ps 

Recall that the coefficient of the pentagon C at one loop was invariant under cyclic 

permutations of the labels. Here the symmetry has been explicitly broken by the extra 

factor of S12 and hence the integral 1^ must be inside the sum over cyclic permutations. 
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FIG. 8: Evaluation of the sum over Feynman diagrams on two different T 8 contours. The residue of 
the p integral on any of the two T 4 's represented by the lines circling the propagators is a tree- level 
four particle amplitude. For n = 5, as opposed to n = 4, there are two choices for the T 4 contour 
in the q variable. Here we choose the one that induces a factorization in the (q — k\) 2 channel of 
the tree amplitude. The other T 4 is also important and is considered in figure [TU1 The residue on 
the final T 4 is the five-particle tree- level amplitude on q^ and zero on . 

2. Second Topology 

The next topology of Feynman diagrams is shown in figure [7(3. The contributing integrals 
are clearly JW and I^ e ' in figure O After performing the integration over p on a T 4 contour 
and choosing the remaining T 4 in the q variables to induce a factorization on the four-particle 
amplitude with 1 and 2, we go back to a one- loop calculation related to the one done in the 
previous case by a cyclic permutation. 

The solution for the coefficients is obtained by performing a cyclic permutation on the 
one-loop coefficients and then multiplying by the I/S12 factor to convert them into two- loop 
coefficients. The explicit form of the coefficients is 

2/5 2 
^2-loop _ _^trce^5l£l2P4 ^2-loop _ ^tree g 45 g 51 g 12 

Pa — 04 Pa — f3 a 
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(A) (B) (C) 

FIG. 9: Scalar integrals contributing to the third kind of topology. In the computation of the 
residues the propagators not used as poles remain and must be evaluated at the location of the 
singularity. There is one propagator left in each pentagon-box integral and they are enclosed by 
dashed lines in the figure. 

Note that we have computed the coefficient of the pentagon-box integral, I^ e \ once again 
and the answer looks very different. Using the same identity (fl2"j) as in the previous section 
it is clear that the two formulas are the same since all we have done is to multiply them by 
the same Si 2 factor. 



3. Third Topology 

Let us consider the final topology shown in figure HC. The possible scalar diagrams that 
contribute to the T 8 which computes the natural leading singularity of this topology are 
shown in[9j Let us denote their coefficients by B, C5 and C3 respectively. The coefficients 
C3 and C5 have already been computed. So in practice one would need a single equation to 
determine B. Let us however consider the two equations that arise and check that they can 
be consistently solved. 

The situation here is different from the situation in the previous cases since the sum over 
Feynman diagrams has zero residue on the leading singularities we will consider. From that 
point of view, we can think of the integral in figure [9j\ as canceling unphysical singularities in 
the integrals in figures 03 and [9p which were already shown to be present in the amplitude. 

The location of the unphysical leading singularities is given by 

„(1) _ [Ml \ \ _ ( 5 > 4 ) \ \ foe \ 

p "M 5 ' q ~ ~ (M) (26) 
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and 

^ = §1^, ,27, 

In order to understand how these arise note that the middle propagator, l/(q — p) 2 , in 
figure [9]A. has two poles on the locus where p 2 = q 2 = 0. These are (p, q) — and [p, g] = 0. 
Using both of them at the same time gives the T 8 used to obtain ( 1261) and ( 1271) . 

Since p and g are proportional to each other in both solutions it is clear that the sum 
over Feynman diagrams shown in figure [7p vanishes. 

The equations we have to solve take the form 

B + (qV-kt-h) 2 + {p^-k 2 -hf = ° (28) 

for i — 1,2. The consistency condition that the known coefficients C% and C5 must satisfy 
for these equations to have a solution is 

9i + 9i = 9i + 9i (29) 

{qW - k x - k 2 ) 2 (pW - k 2 - h) 2 (g( 2 ) -ki- h) 2 (p( 2 ) - k 2 - hf K ' 
Using the explicit form of the coefficients 

c = Atrce £34^51 = Atme S 51 S 2 2 S2 3 

one can check that (129]) is indeed satisfied. Therefore any of the two equations determine 
B. Choosing i = 1 gives 

B = ^5 C 3 

( g (i) _ kl - k 2 ) 2 ( p W-k 2 -k 3 y y ' 

This expression can be dramatically simplified if (129!) is used to solve for C3 in terms of C5. 



The answer turns out to be 



B = -^. (32) 
S12 



This is the coefficient of integral 1^ in figure El 



4- Additional Leading Singularities 

As mentioned in the discussion of the first kind of topology, once the integration over 
the T 4 corresponding to the p variables is carried out, one is left with a sum over Feynman 
diagrams that possesses two four-particle amplitudes (see figure [8]); one which contains 
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particles 1 and 2 and another which contains particles 4 and 5. Previously, we chose the 
remaining T 4 contour such that the amplitude with 1 and 2 factorizes. Now we have to 
check that the leading singularity corresponding to 4 and 5 factorizing also works. 

Here we get four equations corresponding to any combination of and p^\ p^ . 

Note that none of the scalar integrals with the topology of two boxes, i.e. I^ a \ I^ b \ 1^ 
contributes to these leading singularities. At this point, we have a single integral that 
contributes, i.e., I^ e \ Not surprisingly, it is not possible to solve all the equations with the 
coefficient of a single integral. This means that at least one more integral is missing. 

The requirements for the new integral are: 

• It must be non-zero on all the new T 8 contours. 

• It must vanish on all the previous T 8 contours since they have already been accounted 
for. 

The only possibility is to start with an integral of the form 1^ which ensures the first 
requirement and then introduce a zero in the numerator with removes the pole (q — ki) 2 
which enters in the definition of all other T 8 's already studied. Removing the pole at (q — k±) 2 
guarantees that the new integral has zero residue in all previous leading singularities thus 
fulfilling the second requirement. 

The new integral has the topology of 1^ and a numerator factor (q — ki) 2 . This integral 
is precisely 1^ in figure 

Let us show that by adding to the basis all remaining leading singularities can be 
accounted for. Let its coefficient be D. 

The four leading singularities are located at 

p[i){q)= W] XiK p{2){q)= W) Xq ^ (33) 

and 

^-^+114 « m -( x * + m x ¥>- (34) 

The value of p does not enter in the computation and we only have to consider the two 
equations coming from the two values of g's. The equations are obtained in a completely 
analogous way from the previous calculations following the steps in figure [TU1 The two 
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FIG. 10: Evaluation of the sum over Feynman diagrams over the new T 8 contours. In figure [8] we 
chose the second T 4 to induce a factorization in the (q — k±) 2 channel, here we choose it to induce 
a factorization of the second four-particle amplitude in the (q + k$) 2 channel. 



equations are 

CD C 
Si2S 3 4S i5 (q^ - h) 2 + S34S45S12 A ' (<? (2) - h) 2 + D °' ^ 

The coefficient C has already been computed so in practice one would simply use the second 
equation to determine D. Let us not do that and solve the equations for D and C. 
Defining 

, s fi t isr, 7:= fi+MMV (3 6 ) 



(3,5>[5,1]7 ' 'V <5,1)[3,5] 
the solution is given as follows 



c = ^ ^45^1^12 D = _^c Si2S34S45 ^__ (37) 

(7-7) 7-7 

An explicit computation reveals that the new expression for C agrees with that of C 2 ' loop 
found earlier. 
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5. Final Result 



We now collect all the results obtained by requiring that all leading singularities are 
reproduced correctly. Using the scalar integrals defined in figure M the two-loop amplitude 
is 



| (2) 

cyclic 

where 



A 5 = X^ S (_ £12^23^5 j(a) _ S5l£l2^4 _ £34f457 j(c ) _ ^51^23 j( d ) + ^51^12^23 ^) 

4 0) £t\ /?5 -/3 5 /?4-/3 4 7-7 fo-h fc-h 



(38) 



a - = 1 1+ ftmtwmJ ' 7 - = 1 1+ (3^)mJ ' (39) 

and /3j and 7 the parity conjugated expressions of $ and 7 respectively. 

An even simpler expression can be obtained if one uses the different identities found in 
the previous subsections to write 

1 S 12S23 1 1 s 23 1 



7-7 S34S45 (A, -AO' *45 (#>-#>)' 

Using (jlOl in (jggjl one finds 



(40) 



E f-«i2^/ (a) - f^/W - Sl2 7l< c > - S5 i/ {d) + s, lSl2 A . (41) 



In order to compare with the known answer it is convenient to split the coefficients in parity 
even and odd terms. Note that using out technique both pieces are computed simultaneously 
and equally straightforwardly. The separation into definite parity terms is easily done by 
recalling that for real momenta, (3, 7 are the complex conjugate to (3 and 7 respectively. 
This means that we can write, for example, 

- f^/w = f 1 - ^rHr) < 42 > 

P5 — P 5 * \ p 5 — p 5 J 

Note that some coefficients are naturally parity odd. Like those of 1^ and /w. 

It is easy to check analytically using a symbolic manipulation program like Mathematica 
that, up to an overall normalization of 1/4, ( 1381) is exactly equal to ( 1201) . 
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VI. A PEEK AT TWO-LOOP SIX-PARTICLE AMPLITUDES: MHV AND NEXT- 
TO-MHV 

One of the advantages of our technique is that the homogenous part of the linear equations 
that determine the coefficients of the integrals is helicity independent! All the helicity 
information enters in the inhomogeneous part. 

In order to illustrate this feature the first non-trivial case is that of six particles. This 
is the first case where next-to-MHV configurations are possible. In this section we choose 
a particular subset of Feynman diagrams contributing to two-loop six-particle amplitudes. 
Studying the leading singularities one can write down linear equations that determine the 
complete coefficient of all pentagon-pentagon integrals and of a certain class of pentagon- 
box integrals. The determination of the complete set of linear equations which gives the full 
amplitude is outside the scope of this paper and we leave it for future work. 

The set of Feynman diagrams we consider generates the five topologies shown in figure 
rm We write down explicitly the equations coming from (A) and (B), as the remaining 
three, (C), (D), and (E), are completely analogous to (B). 



A. Topology (A) 

The integrals contributing to the first kind of topology are shown in figure O If we were 
to follow the same steps as in the previous section we would start by taking only the first two 
integrals in the figure and then realize that it is not possible to solve the four equations that 
arise by comparing all leading singularities. From the experience with the five-particle case, 
we start by adding eight integrals with numerators such that they will contribute to the four 
leading singularities under consideration but they will not contribute to other singularities 
where they are not needed just like in the case of j' c ) for five particles. 

Using the labels in the figure, the four leading singularities are found by choosing any 
combination of p*'s and g*'s from 

(i) _ (2, 3) ~ ( 2 ) _ [2, 3] 7 (i) _ (5, 6) - (2) _ [5, 6] - 

v -<i^y A ^ p ~]ij] x ^ q --Ja^) X ^ q -"M 5 { ) 

Let us denote the evaluation of the sum over Feynman diagrams (see figure [TTK) in a 
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FIG. 11: Topologies of sums over Feynman diagrams which can be used to determine the coeffi- 
cients of all pentagon-pentagon and some class of pentagon-box integrals in a six-particle two-loop 
amplitude. Some choice of external labels has been made. No helicities have been assigned since 
all choices, MHV and next-to-MHV, can be treated simultaneously. 

particular solution (p, q) by F Ptq . Then the equations are 



P.'l 



B + 



{p+q \ 2 ^ (C + (£1 + D 1 (q - k 2U ) 2 )(p + k 61 ) 2 + 



(E 3 + D 2 (p - A; 234 ) 2 )(g - k 23A ) 2 + (E 4 + D 3 (p + A; 61 ) 2 )(g - k u ) 2 + 
(E 2 + D 4 (q-k u ) 2 )(p-k 2U ) 2 )] 

(44) 

The labeling of the coefficients is explained in the caption of figure [121 

The system at hand involves four equations and ten unknown coefficients. In order to 

find a set of equations sufficient to completely determine the coefficients we have to consider 

the other topologies in figure dU 

Before going to the next topology, lead us compute F(p, q) in some cases in order to 

illustrate the procedure which turns out to be fairly simple since the computation is reduced 

to that at one-loop. 
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FIG. 12: Scalar and generalized scalar integrals that can contribute to the first kind of topology. 
Their coefficients are B, C, Ei, E2, E3, E4, D\, D2, D3, D4 respectively. In order to avoid cluttering 
of the figure we have added the minimal amount of information needed to determine the labels. In 
the double box diagram, other external legs are labeled following the color ordering. All pentagon- 
pentagon diagram naturally inherit their labeling from the double-box diagram. 

All the steps are shown in figure [T3J Consider the lower box in the two-loop diagram. 
These are the Feynman diagrams in a one-loop five-particle amplitude where two of the 
external momenta are actually internal legs in the two-loop diagram. From the one-loop 
discussion in the previous section we know that depending on the helicities, one solution, 
p*, gives zero while the other gives A^ ee . Plugging this into the two- loop diagram one gets, 
in the non-zero case, a one-loop six particle diagram. If the amplitude we are trying to 
compute is MHV or MHV then the answer is either zero or v4g ree . 

As an explicit example take helicities to be {1~, 2~, 3 + , 4 + , 5 + , 6 + }. Then we find 

F(p ( - 1 \ q^) = A^ e , F(pW, q {2) ) = F{ P M, q^) = F(p®, q {2) ) = 0. (45) 



27 



next-to-MHV (split) 







FIG. 13: Computation of functions F Pi g in two cases: MHV with helicity configuration 
{l-,2-,3 + ,4+,5 + ,6 + } and next-to-MHV (split), i.e. {1~, 2~, 3 - , 4+, 5+, 6+} 

If the amplitude is next-to-MHV then the answer is more interesting as it corresponds 
to the coefficient of a particular one-mass integral in a six-particle one-loop amplitude. All 
six-particle one- loop next-to-MHV amplitudes amplitudes where computed in the early 90 's 



in 



They are all given in terms of a quantity defined in eq. 6.13 of [8J] (we have parity 



conjugated the expression in [8j) 



(l|2 + 3|4](3|l + 2|6]tf 23 



°" [l,2][2,3](4,5)(5,6)(ti23t345 -Si2S45)(tl23*234 -S23S56)' 

As a explicit example take helicities to be {1~, 2~, 3~, 4 + , 5 + , 6 + }. Then we find 

F(p« g (1) ) = F(p®, g«) = F{jP\ qW) = 0, F(pU, q^) = B . (47) 
It should be clear that any other helicity configurations can be treated analogously. 



B. Topology (B) 

The scalar integrals that contribute to these leading singularities are shown in figure ITU 
The location of the leading singularities is slightly more involved as p is determined as a 

function of q. Let us give the two solutions for q and then give the two g-dependent solutions 

for p. 
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FIG. 14: Scalar and generalized scalar integrals that can contribute to the second kind of topology. 
Their corresponding coefficients are F, G, C, E\, E2, E3, £4, D\, D2, D$, D4 respectively. Just as in 
figure [T2l the labels in the diagram are determine from the those of the first integral. 



and a,P the parity conjugate of a, (3. 

Once again we can easily write down the linear equations coming from imposing the 
correct behavior at the four leading singularities. Let us denote by H p q the value of the sum 
over Feynman diagrams on one of the four solutions p*, q*. Then the equations are 



Wi2S 23 (<? - h6) 2 H m = F + G(q - k 56 ) 2 + ? --Ur (C + (E 1 + D 1 {p - h) 2 )(q - h) 2 + 



The labels of the coefficients are explained in the caption of figure [TH 

Combining the two systems of equations, (jH| and (15D|1 . one finds eight equations. 
Repeating the same procedure for the three remaining topologies in figure [TT] which are 

completely analogous to topology (B) one finds four equations in each case. This gives a 

total of twenty equations for seventeen coefficients, i.e., eight pentagon-box integrals and 

nine pentagon-pentagon integrals. 
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Using seventeen equations to determine the coefficients leaves three more as consistency 
checks. 

We end this section by mentioning that the next natural set of equations to consider 
comes from double-box topologies where one of the boxes is only a four-particle one-loop 
amplitude. For a given choice of the two external legs to the four-particle amplitude, say 
{1,2}, there are three topologies. These correspond to the possible ways of distributing 
the remaining external legs; in our example we would find {3, 4, 5; 6}, {3, 4; 5, 6}, {3; 4, 5, 6}. 
By looking at the leading singularity which uses the propagator coming form the jacobian 
of the four-particle box, these give rise to only two equations each. Hence we find six 
equations in total. There are however seven scalar integrals that contribute. One of them 
is a hexagon-box. 

Following this line of ideas it would be interesting to continue and write down the complete 
set of linear equations which determines all two-loop six particle amplitudes. 

Already with the results presented here, it would be interesting to compare the parity 
even part of the coefficient of all pentagon-pentagon integrals to those obtained recently 
in 27|] for MHV amplitudes. 



VII. CONCLUSIONS 



The program of determining the S-matrix of physical theories by studying the structure 
of its singularities when analytically continued into complex values of kinematical invariants 
was very ambitious. Now we know that one of the main goals of the program which was to 
understand the strong interactions has proven to be a formidable problem. 

In the past 20 years we have learned that M = 4 super Yang-Mills (SYM) can serve as a 
laboratory where new techniques and ideas can be tested. Based on the results presented in 
this paper we would like to propose that perhaps M = 4 SYM, at least in the large N limit, 
can be used to realize the basic ideas of the S-matrix program. Of course, M = 4 SYM 
in four dimensions is a conformal theory and no ^-matrix can be defined. However, if IR 
divergencies are regulated using dimensional regularization then a sensible 5*-matrix can be 
defined. It is very intriguing that on the contours which computer the discontinuity across 
leading singularities, all integrals are finite and hence make perfect sense in four dimensions. 
It is natural to expect the amplitudes on the torus contours can have physical meaning. It 
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would be interesting to explore this further. 

In general, Feynman diagrams possess many singularities; from poles, related to reso- 
nances, to branch cuts, related to unitarity. When a given scattering amplitude in theories 
with spin is considered as a function of independent complex variables Xa^ and the 
problem of reconstructing it from the structure of its singularities is certainly out of hand. 



At tree-level we learned few years ago 
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351 ] that one can solve a simpler problem by 



concentrating on a single complex variable deformation of the amplitude. The main sim- 
plification comes from the fact by matching only a very small subset of all poles is enough 
to determining the amplitude in terms of smaller ones. This idea led to recursion relations 
between on-shell scattering amplitudes. 

At loop level, Feynman diagrams posses an intricate structure of nested branch cuts. It 
turns out that the discontinuity across a given branch cut possesses branch cuts itself. The 
nested structure gets more an more complicated the higher the loop order. The problem 
of finding functions which reproduce all such discontinuities is clearly very hard. Out of 
all these singularities, the special class studied in this paper, called leading singularities, 
are the ones which have the highest codimension. We have argued that the discontinuities 
associated to them are determined by simple computations of residues. We have shown that 
quite remarkably, the problem of finding a function which reproduces all leading singularities, 
which only implies the solution to linear equations, is enough to completely determine the 
amplitude in J\f = 4 SYM. 

Summarizing, computing multi-loop, multi-particle amplitudes in M = 4 SYM can be 
reduced to the computation of residues and the solution of systems of linear equations. The 
residues only involve products of tree-level amplitudes. These possess all the helicity infor- 
mation and only enter in the inhomogeneous part of the linear equations. The homogeneous 
part is universal. If we think about the amplitude and the linear equations as being ob- 
jects with the same information, then it might be that the matrix which determines the 
homogenous part of the equations is the right object to study properties like strong coupling 
expansions, collinear limits, IR consistency equations, etc. 

A natural question that arises is whether this is property is unique to M = 4 SYM. The 
most promising theories are those for which one-loop amplitudes can be written in terms of 
only boxes, i.e., bubbles and triangles are absent. Af = 8 supergravity has been hypothesized 



to have such property 
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3, 



201 ] . This means that one could try to apply the technique 
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presented here to that case. Already, an important step in this direction was given in 16]. 
It would be very interesting to compute more examples. 

In 16J], higher loop four-particle amplitudes in TV = 4 SYM were considered. In cases 



where four-particle one-loop amplitudes were present as subdiagrams it was shown that 
one could related the computation to that at one less loop order. It would be interesting 
to explore the consequences that matching all leading singularities would impose on the 
generalized scalar integrals. Perhaps a formal derivation of the idea of corrections introduced 
in [3] can be found. 

As a final note one can say that the power of matching each individual leading singularities 
comes from the fact that the number of conditions which are naturally linear grows much 
faster with the number of loops than those of previous approaches. 
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